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We discuss the octupole deformation properties of the recently proposed Barcelona-Catania-Paris 
(BCP) energy density functionals for two set of isotopes, those of radium and barium, where it is 
believed that octupole deformation plays a role in the description of the ground state. The analysis 
is carried out in the mean field framework (Hartree- Fock- Bogoliubov approximation) by using the 
axially symmetric octupole moment as a constraint. The main ingredients entering the octupole 
collective Hamiltonian are evaluated and the lowest lying octupole eigenstates are obtained. In this 
way we restore, in an approximate way, the parity symmetry spontaneously broken by the mean 
field and also incorporate octupole fluctuations around the ground state solution. For each isotope 
the energy of the lowest lying l~state and the B(E1) and B(E3) transition probabilities have 
been computed and compared to both the experimental data and the results obtained in the same 
framework with the Gogny D1S interaction, which are used here as a well established benchmark. 
Finally, the octupolarity of the configurations involved in the way down to fission of 240 Pu, which is 
strongly connected to the asymmetric fragment mass distribution, is studied. We confirm with this 
thorough study the suitability of the BCP functionals to describe octupole related phenomena. 
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I. INTRODUCTION 

The ground state and low lying excited states of many 
atomic nuclei all over the nuclide chart show quadrupolc 
deformation in their intrinsic states This property 
has profound consequences in the low lying spectrum 
of those nuclei as well as in their decay patterns [U, HJ ■ 
Octupole deformation is not as common as quadrupole 
deformation as a characteristic of the ground state of 
atomic nuclei, but its consequences are important for un- 
derstanding nuclear properties of several actinide nuclei 
around radium and several rare earth around barium [J . 
The octupole operator has negative parity and therefore 
a non-zero octupole deformation means that the intrinsic 
state has lost reflection symmetry and acquired a pear- 
like shape. The quantum interference between the two 
degenerate intrinsic states with pear-shaped matter dis- 
tributions pointing upward and downward (i.e. with the 
same absolute value of the octupole moment but oppo- 
site sign) restores the parity quantum number and leads 
to the presence in the spectrum of a doublet with op- 
posite parities [U, 0]. The energy splitting between the 
two members of the doublet strongly depends upon the 
properties of the barrier separating the two degenerate 
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intrinsic states with opposite octupole moment. In de- 
formed even-even nuclei it is possible that the negative 
parity member of the ground state multiplet , the lowest 
lying 1~ state, can be located below the lowest 2 + state 
leading to the appearance of alternating parity rotational 
bands, which are clear signatures of octupole deforma- 
tion. Also, the two members of the doublet will be con- 
nected by strong B(E1, 1~ — > + ) transition probabili- 
ties from the 1 _ to the ground state. The next mem- 
ber of the negative parity rotational band is a 3~ that 
rapidly decays to the + ground state by means of strong 
B(E3, 3~ — > + ) transition probabilities. Although there 
are several known examples of alternating parity rota- 
tional bands at low spins, the alternating behavior usu- 
ally appears at high spins as a consequence of the sta- 
bilizing effect of angular momentum on the octupolarity 
of the system: moments of inertia increase with the oc- 
tupole moment and therefore configurations with higher 
octupole moments are the more lowered by increasing 
angular momentum. 

The appearance of octupole effects is strongly linked 
to the position of the Fermi energy in the single parti- 
cle spectrum of the underlying mean field [l|, The 
reason is that octupolarity is enhanced when, in a given 
major shell, the intruder orbital interacts (via a particle- 
hole excitation) with a nearby normal parity orbital with 
three units less of angular momentum (3h is the amount 
of angular momentum carried out by the octupole oper- 
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ator). This happens between the j'15/2 and g 9 / 2 spheri- 
cal orbitals, the I13/2 and f 7 / 2 and the h-^n and d 5 / 2 - 
Those regions where both protons and neutrons feel a 
strong octupole interaction is where octupole related ef- 
fects are expected to be more pronounced. For example, 
in the region around |g Ra, the Fermi level of protons 
is located around the /^orbital that can then interact 
strongly with the empty 113/2- Besides, Fermi level of 
neutrons is around the g g / 2 orbital, which strongly inter- 
acts through the octupole interaction with the j 15 / 2 one. 
Similar arguments apply to the region around 144 Ba with 
the orbitals d$/2 and /in/a for protons and the /V/2 and 
H3/2 f° r neutrons responsible for octupolar effects. 

Last but not least, octupolarity also plays a relevant 
role in the asymmetric fission decay mode because an 
octupole deformed pathway to fission naturally explain 
the observed asymmetric mass fragment distribution of 
several actinide parent nuclei [5|. 

In this article we want to check the ability of the re- 
cently proposed Barcelona-Catania-Paris (BCP) energy 
density functionals [6], in dealing with the octupole de- 
gree of freedom in finite nuclei. We will compare the re- 
sults provided by the BCP functionals in some test cases 
with experimental values, when available, and with the 
results obtained using the Gogny D1S interaction that 
we take here as a benchmark. The BCP energy density 
functionals consist of a bulk part which is fully micro- 
scopic and comes from the nuclear and neutron equa- 
tion of state which are parametrized in a polynomial 
form complemented by additional terms accounting for 
finite size effects (see [3, |9j for functionals inspired in 
the same principles). In addition to the Coulomb term 
and the spin-orbit contribution, which is taken exactly 
as in the Skyrme or Gogny forces, we add a purely phe- 
nomenological finite range term for describing properly 
the nuclear surface. To deal with open-shell nuclei we 
still include in the BCP functionals a zero range density- 
dependent pairing interaction fitted to reproduce the nu- 
clear matter gaps obtained with the Gogny force (l2^ . 
The only free parameters of these functionals are the 
isospin like (L) and unlike (U) strengths of the surface 
term, the range of the Gaussian form factor used to give 
a finite range to the surface term and the strength of the 
spin-orbit interaction Q. These free parameters are ad- 
justed in the usual way to reproduce the ground state 
energy and charge radii of some selected spherical nu- 
clei. With these ingredients, the BCP functionals give 
an excellent description of 161 even-even spherical nuclei 
with mis values for the ground-state energies (1.77 MeV 
and 2.06 MeV) for BCP1 and BCP2, respectively) and 
charge radii that are comparable to the ones obtained 
with well reputed interactions/functionals like Skyrme 
SLy4 (1.71 MeV), Gogny D1S (2.41 MeV) or the Rel- 
ativistic NL3 parametrization (3.58 MeV). Apart from 
the advantages already mentioned in the BCP func- 
tionals are advantageous in its application to finite nuclei 
because of its reduced computational cost as compared 
to Gogny (BCP is a factor between 6 to 10 faster) or 



even Skyrme (comparable computational cost). Also the 
appearance of integer powers of the density in the bulk 
part of the functional, which is a consequence of the spe- 
cific fit to the nuclear matter results, makes much easier 
to deal with the self- energy problem that plagues beyond 
mean field calculation [10| . Using these BCP functionals 
we have also explored quadrupole deformation proper- 
ties [ll|. We find a behavior similar to that obtained 
using the Gogny D1S force widely used to this end. This 
fact give us confidence in using the BCP functionals to 
study nuclear properties related to deformation. As the 
BCP functionals are aimed at describing not only masses 
and radii but also the low lying spectrum over all the 
nuclide chart, it is necessary to check whether the very 
reasonable results regarding quadrupole collectivity can 
also be extended to the octupole deformation case. To 
check that this is the case, we have carried out mean field 
Hartree- Fock- Bogoliubov (HFB) calculations with the 
BCP energy density functional as well as the Gogny [l3| 
D1S [HI interaction to test the response of the system to 
the octupole degree of freedom. To be more precise, we 
have used a constraint in the axially symmetric octupole 
moment to generate potential energy curves (PECs) to 
search for octupole deformed minima as well as to study 
the stiffness of those (and other) minima against changes 
in the octupole degree of freedom. These PECs are com- 
puted for several isotopes of radium from 216 Ra to 232 Ra 
and of barium from 140 Ba until 150 Ba. In addition to the 
PEC, the calculation of the corresponding collective iner- 
tias allows the evaluation of the 1~ excitation energy as 
well as B(E1) and B(E3) transition probabilities in the 
framework of the collective Schrodinger equation (CSE) 
method. The results will be compared with experimental 
data, when available, as well as with the results obtained 
with the Gogny D1S force. It should be mentioned that 
the Gogny D1S results have already been reported in 
Refs. [15l4l7| and similar calculations with the Skyrme 
interactions exist in the literature [18j. Finally, the oc- 
tupole properties of the fission valley of 240 Pu will also be 
briefly discussed and compared to those of Gogny D1S. 



II. THEORETICAL TOOLS 

To solve the Hartree-Fock-Bogoliubov (HFB) equation 
the quasi-particle operators of the Bogoliubov trans- 
formation have been expanded in a harmonic oscillator 
(HO) basis big enough as to warrant convergence of the 
results with the basis size. The expansion coefficients 
have been determined by means of the gradient method 
which relies on the parametrization of the mean field 
(HFB) energy in terms of the parameters of the Thou- 
less expansion of the most general HFB wave functions. 
Within the gradient method, the HFB problem is recast 
in terms of a minimization (variational) process of the 
mean field energy and the search of the minimum is per- 
formed by following the direction of the gradient in the 
multidimensional space of parameters. The advantage 
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of this method over the more traditional one of succes- 
sive diagonalizations is in the way the constraints are 
implemented, which allows a larger number of them to 
be treated at once. Axial symmetry has been preserved 
in the calculation implying the use of an axially symmet- 
ric HO basis made up of the tensor product of two di- 
mensional HO wave functions times one dimensional HO 
ones. Along with the octupole moment constraint asso- 
ciated to the multipole operator Q3 = r 3 l3o and used 
to generate the PEC's, we have included a constraint on 
the center of mass of the nucleus (i.e. the mean value 
of r 1 Yio has been set to zero) to prevent spuriousness 
associated to the center of mass motion, to slip into the 
results. As a consequence of the axial symmetry imposed 
in the HFB wave functions, the mean values of the 
multipole operators with fi ^ are zero by construction. 

The information given by mean field theories is re- 
stricted to the energy and shape of the -generally- de- 
formed ground state. To restore the parity symme- 
try broken by the mean field approximation and to de- 
scribe the dynamics of the collective excited states it 
is mandatory to go beyond the mean field approxima- 
tion. With this in mind, the octupole degree of freedom 
Q3 = {ip\Q3\ip) (where \ip) is the HFB intrinsic wave 
function) has been used to build up a collective Hamilto- 
nian based on the Generator Coordinate Method (GCM) 
and the Gaussian Overlap Approximation (GOA) [19l — 
l2ll |. In this method, the GOA is used to reduce the Hill- 
Wheeler equation of the GCM to a Schrodinger equation 
for the collective wave function, the so-called Collective 
Schrodinger Equation (CSE) 

/ Hcoll4'a(Q3) — ta4>a{Q?,), (1) 

where the collective Hamiltonian TA. C oU is given by 



+ V(Q 3 )-e Q (Q 3 ). (3) 

In this expression G(Qs) is the metric, B(Qz) is the 
mass parameter associated with the collective motion 
along Q3, V{Qz) is the collective potential given by the 
HFB energy V(Q 3 ) = (fl>(Qs)\A\ip(Q 3 )) and e (Q 3 ) is 
the Zero Point Energy (ZPE) correction. The eigenfunc- 
tions <t> a {Qz) of Eq. ([1} have to be normalized to one 
with the metric G{Q 3 ) 

j 'dQ 3y /G(Q3)<l>* a (Q 3 )MQ3) = 6 a ,P ( 4 ) 

to preserve the hermiticity of H C oii- 

It should be mentioned that a CSE can also be ob- 
tained from the Adiabatic Time Dependent Hartree- Fock 
(ATDHF) theory [H-SH after quantization of the semi- 
classical Hamiltonian for the slow moving collective de- 
grees of freedom. The collective Hamiltonian obtained in 
this way has the same functional form as the GCM+GOA 
one, but the expression of the collective parameters is 



different. Later we will discuss how to choose these col- 
lective parameters. 

An interesting characteristic of the collective Hamilto- 
nian for the octupole degree of freedom is that H C oii is in- 
variant under the exchange Q3 — > —Q3 and, therefore, it 
is possible to classify its eigenf unctions, (j> a {Q 3 ), accord- 
ing to their parity under the Q3 — ¥ —Q3 exchange. It is 
easy to see that the parity of the collective wave function 
under the Q3 — > —Q3 exchange corresponds to the spa- 
tial parity operation in the correlated wave function built 
up from 4> a . The inclusion of octupole correlations im- 
mediately restores the parity symmetry lost at the mean 
field level. Therefore, the solution of the CSE Eq. flTJ al- 
lows the calculation of the + — 1 _ (3~) energy splitting 
and the B(E1) and B(E3) transition probabilities con- 
necting them. At this point it has to be pointed out that 
in the present framework where only time reversal invari- 
ant wave functions are considered it is only possible to 
describe excited states with average angular momentum 
zero. To deal with genuine 1 _ or 3 _ states cranking model 
wave functions should be considered, which is out of the 
scope of the present work. Here we will assume that 
the cranking rotational energy of the 1~ state is much 
smaller than the excitation energy of the negative parity 
band head and therefore can be safely neglected. Also 
the impact of the cranking term in the transition prob- 
abilities to be discussed next is neglected. With these 
approximation in mind, the reduced transition probabil- 
ities from the lowest l - and 3~ states to the + ground 
state can be computed within the Rotational Model ap- 
proximation as 

B(EX, If -> /;) = e 2 (/ i ifAO|7 / ir) 2 |( W |r A y A ,o|^/)| 2 , 

(5) 

where \ipt) and \y>f) are correlated wave functions ob- 
tained in the spirit of the GCM from the collective wave 
functions <f> a {Qs,)- The preceding formula can be reduced 
to an expression involving those collective wave functions 
0a (Q3) by means of the GOA (25[. The final result for 
K = bands reads 

B(El, 1- 0+) = |- |(0o- |A)|0 o +)coll| 2 (6) 
for the El electric transition and 

B(E3,3- -> 0+) = |- |(0o- IQ3o(PROT)|0 o+ )coll| 2 , 

(7) 

for the E3 one. In the preceding formulas we have in- 
troduced the collective matrix element of an operator O 
as 

(00- |O|0 O +)COLL = / dQsG 1 / 2 0S-(Q 3 )O(Q 3 )0o+(Q3) 

where 0(Q 3 ) = (^(Q 3 )|6|^(Q 3 )). In Eq. © D is the 
dipole moment operator whose mean value is defined as 
the difference between the center of mass of protons and 
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neutrons 



III. RESULTS 



N Z 
D (Q 3 ) = ^(Qa^otMQsV-jMQsy^MQs)). 

(8) 

Finally, Q 3 o (PROT) is the part of the octupole operator 
acting on proton's space. 

To carry out the collective calculations it is necessary 
to specify the collective parameters G(Q 3 ), B(Q 3 ) and 
ea{Q^) appearing in the definition of T-Lcoii Eq. @. As it 
was said before, there are two sets of parameters coming 
from the GCM+GOA and the ATDHF derivation of the 
collective Hamiltonian. The set of parameters used in 
this calculation is an admixture of the two and it is known 
as the ATDHF +ZPE set. It includes the mass parameter 
B(Q 3 ) coming out from the semi-classical Hamiltonian 
of the ATDHF theory, the metric of the GCM+GOA 
and the ZPE correction computed with the GCM+GOA 
formula but using the ATDHF mass instead, i.e. 
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G(Q 3 )B(Q 3 ) 



ATDHF- 



(9) 



This set of parameters was devised to put together the 
advantages of the ATDHF set (time-odd components in- 
cluded in the mass term) and the ones of the GCM+GOA 
(ZPE correction). This method can be somewhat justi- 
fied in the context of the extended GCM |2ll . |26| and has 
been extensively used [H, [l6[ . 

The calculation of the collective parameters involves 
the inversion of the HFB stability matrix which is closely 
related to the matrix of the RPA equation. At present, 
this is a formidable task and approximations are needed. 
The approximation used in this article - called "cranking 
approximation" [27ll28j - neglects the off-diagonal terms 
of the stability matrix allowing to invert it analytically 
but at the cost of including the two body interaction 
only through the mean field. Although this approxima- 
tion has been extensively used in the literature for the 
calculation of collective masses and moments of inertia 
(see, for instance, Refs. [3, H^, [3(j ) its validity has not 
been properly established. Using the "cranking approxi- 
mation", the ATDHF+ZPE parameters are given by 



G(Q 3 ) 



M- 2 (Q 3 ) 
2M 2 _ 1 (Q 3 y 



B(Q 3 ) 



M- 3 (Q 3 ) 
M\{Q 3 Y 



(10) 



where the quantities M_ n (Q 3 ) (n — 1,2,3 ) are defined 
as 



M_„(Q 3 ) = 



\(Q30)t°i 



k.l 



(11) 



In the preceding expression, Ek are the quasi-particle 
energies and {Q 3 o) kl are the matrix elements of the 20 
part 0| of the octupole operator Q30 in the quasi-particle 
basis of the HFB wave function \ip(Q 3 )). This form of the 
collective mass is usually referred to in the literature as 
Belyaev-Inglis mass [2] . 



In the subsequent subsections the results obtained with 
the BCP1 functional and regarding octupole prop- 
erties of some radium and barium isotopes will be dis- 
cussed. The other functional defined in [6] and referred to 
as BCP2 will not be explicitly considered here although 
the calculations were carried out for that case too. The 
reason is the strong similarities between BCP1 and BCP2 
results that produced most of the curves one on top of 
another, making it impossible to differenciate them in the 
plots presented. 



A. Low excitation energy properties in the radium 
isotopes 

Octupole deformation properties of the radium iso- 
topes were the first to be addressed from a microscopic 
point of view with the Gogny force, first at the mean 
field level [l6[ and next including exact restoration of the 
parity symmetry [3l|. At the mean field level, the first 
quantity to analyze is the potential energy curve (PEC) 
as a function of the octupole moment that determines 
both the ground state minimum and its stiffness. Let 
us point out that for every point in the PEC the other 
multipole moments (quadrupole, hexadecapole, etc) are 
selfconsistently determined as to produce the lowest en- 
ergy. The PECs computed with the BCP1 [f| energy 
density functional and the Gogny D1S [l3|, [l4| force are 
depicted in Fig. [T] As can be seen in the plot, the re- 
sults for the two types of interactions look very similar 
in all the nuclei considered. It is observed how, whenever 
a minimum appears (in the nuclei from 218 Ra to 228 Ra) 
in the Gogny D1S calculation at an octupole deforma- 
tion different from zero, the same happens and at the 
same Q 3 value in the BCP1 calculation. For the nuclei 
with the minimum at Q 3 = the Gogny D1S force shows 
a tendency to produce a stiffer parabolic behavior in the 
PEC than in BCP1. The depth of the octupole deformed 
minima is also very similar for both kinds of calculations 
and reaches its maximum value of 1.5 MeV for the nucleus 
222 Ra which therefore can be considered as the strongest 
octupole deformed nuclei of the considered chain. 

In Fig [5] we show the particle-particle correlation en- 
ergy defined as E pp — ^TiAn and given in terms of the 
usual pairing field A and pairing tensor k of the HFB 
method. This quantity gives a rough idea of the amount 
of pairing correlations in the system. It can also be used 
as an indicator of the size of the single particle level den- 
sity around the Fermi surface as strong pairing correla- 
tions are direct consequence of high level densities. This 
energy is also correlated with the pairing gap that rep- 
resents the energy of the lowest two quasi-particle exci- 
tation and therefore it is closely related to the collective 
inertias to be discussed below. The overall tendency of 
Epp is to decrease with increasing octupole moment up 
to values of Q 3 = 10 b 3 / 2 which correspond to typical 
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Figure 1: The HFB mean field energy as a function of the 
octupole moment Q3 (in units of 6 3 ^ 2 = 10 3 fm 3 ) for the 
isotopes of radium (Z=88) from A=216 up to A=232. Results 
for both the BCP1 energy density functional (full curves) and 
the Gogny D1S force (dotted curves) are shown. 




Figure 2: The particle-particle (pp) correlation energy de- 
fined as E pp = ^TrAfi; is plotted as a function of the octupole 
moment Qz (in units of b 3 ^ 2 = 10 3 fm 3 ) for the radium iso- 
topes considered. Dotted lines correspond to the pp corre- 
lation energies whereas full lines are meant to represent the 
same quantity for protons. Thick lines are used to depict the 
results of the calculations with the BCP1 functional whereas 
thin lines are for the results with the Gogny D1S force. 



excitation energies of 5-6 MeV above the ground state 
in the PECs. From there on we observe, depending on 
the nucleus, stationary behaviors or mild increases. We 
also notice that the E pp computed with BCP1 are greater 
than the ones computed with D1S for the light isotopes 
216 Ra and 218 Ra and for the two species of nucleons. For 
the nucleus 220 Ra the particle-particle correlation ener- 
gies for protons and neutrons are similar in both calcu- 
lations and from there on and up to the 230 Ra isotope 
the D1S correlation energies are larger than the BCP1 
ones. For the heaviest isotope considered 232 Ra the E pp 
energy for neutrons is larger for BCP1 than for D1S and 
the opposite is true for protons. 

In Fig. [3] the single particle energies (SPE) are plotted 
as a function of the octupole moment Q3 for both protons 
and neutrons. Results of both calculations are shown in 
different panels ( Gogny D1S upper and BCP1 lower pan- 
els). In all the cases, the Fermi level is represented by a 
thick dotted line. Once the stretching effect of the big- 
ger effective mass of BCP1 (1, versus 0.7 in Gogny D1S) 
is accounted for, the similitude between the single parti- 
cle energies obtained with Gogny D1S and BCP1 around 
the Fermi level and Q3 values in the neighborhood of the 
mean field minimum is remarkable. It is also possible to 
say that the behavior of most of the levels as a function 
of both Q2 and Q3 is quite similar in both kinds of calcu- 
lations. This is so in spite of the different ordering of the 
spherical orbitals: for protons the fa/2 and hg/ 2 are re- 
versed in the spectrum of Gogny D1S as compared to the 
spectrum of BCP1. For neutrons there is also such an in- 
version between the j'15/2 orbital and the in/2 orbital and 
also the separation between the in/2 & ncl 59/2 orbitals is 
much larger in BCP1 than in Gogny D1S. That the sin- 
gle particle spectrum looks similar in the region of well 
developed quadrupole deformation and also as a function 
of octupole deformation is probably a consequence of the 
collective character of those collective degrees of freedom 
where the geometry of the shape of the nucleus is more 
important than quantum mechanical effects. To make 
the argument more quantitative, we have analyzed the 
structure of the single particle wave functions in terms 
of Nilsson quantum numbers and found that the levels 
around the Fermi surface have similar structures. A typ- 
ical example for protons is the K — 5/2" level that for 
Q 3 = lies at -4.5 MeV in the D1S case and at -4.2 MeV 
in the BCP one. This level originates in the h 9 / 2 spheri- 
cal orbital in the D1S case whereas it comes from a fa/ 2 
in the BCP one. The Nilsson quantum numbers for the 
D1S orbital are [523] (66%), [532] (11%), [503] (7%), [743] 
(6%), and other small components whereas for BCP1 
they are [523] (46%), [532] (15%), [503] (9%), [312] (7%) 
and smaller components. We can also consider another 
example in the neutron side where, at Q3 = there is a 
K = 3/2 + orbital at around -5 MeV that originates from 
a ii3/2 spherical level in the D1S calculations and from 
a d 5 /2 in the BCP one. The Nilsson quantum numbers 
obtained are [631] (31%), [642] (25%), [611] (15%), [862] 
(8%) and small components for D1S and [642] (47%), 
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Figure 3: (Color online) Single particle energies (SPE) of the nucleus 224 Ra both for protons (left panels) and neutrons (right 
panels) and the two interactions/functionals used in the calculations (Gogny D1S in the upper panels, BCP1 in the lower 
ones). In each panel, the single particle energies are plotted first as a function of the mass quadrupole moment Q2 starting at 
sphericity (Q2 = b) up to the value corresponding to the self-consistent minimum in the quadrupole variable at Q3 = b 3//2 
(Q2=8.12 b for the Gogny calculation and 7.88 b for the BCP1 one). In this plot the spherical quantum numbers obtained at 
Q2 — b are given. On the right hand side of each panel, the SPE as a function of the octupole moment starting at Q- A — 

b 3/2 

are plotted. The octupole deformed self-consistent minimum is locate in this nucleus at around Q3 = 4 b 3//2 . Thick dotted 
lines indicate the Fermi level. As a consequence of axial symmetry, the K quantum number of each level remains fixed along 
the plot. Each K value has an associated color (black, K = 1/2; red, K = 3/2 ; green, K = 5/2; blue, K = 7/2; dark green, 
K — 9/2, etc). The parity of each state when Q3 = (i.e., in the Q2 plot) is determined by the kind of line: full for positive 
parity and dashed for negative parity. Finally, the K quantum numbers of specific levels around the Fermi level are given. 



[631] (16%), [862] (13%), [422] (7%) and small compo- 
nents for BCP. From the preceding examples and other 
orbitals considered (but not displayed here) we conclude 
that the quantitative structure of the levels is quite sim- 
ilar in the two calculations irrespective of their spherical 
origin. This reinforces our suggestion about the funda- 
mental role played by the collective degrees of freedom 
in the determination of single particle wave functions. 

The conditions for the development of octupolarity are 
clearly satisfied in this nucleus as can be seen in the sin- 
gle particle plot: for protons there are "jy/2" levels below 
the Fermi level with K = 1/2 and 3/2 and at the same 



time orbitals with K = 1/2 and 3/2 are just above 
the Fermi level. The same happens in the neutron case 
with the g$/2 orbital well below the Fermi level and the 
J15/2 with K = 1/2 and 3/2 at the Fermi level (please re- 
member the super-fluid character of neutrons that makes 
the Fermi level concept a diffuse one) . Another condition 
for the development of a minimum is the presence of a 
region of low level density in the SPE spectrum (Jahn- 
Teller effect, see Ref. [1] for a general discussion in the 
nuclear context). We observe in the two proton spectra 
in Fig. [3]how the Fermi level of protons lies in the middle 
of a low level density region at Q3 = 4 b 3 / 2 which corre- 
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sponds to the position of the minimum. For neutrons and 
around Q3 = 4b 3 / 2 we also observe a region of low level 
density near the Fermi level which is more pronounced 
for the BCP1 results. As the number of neutrons is in- 
creased the Fermi level moves upward and enters a region 
of high level density that is unable to lead to a deformed 
minimum as is the case for 230 Ra and heavier isotopes. 

Finally, we would like to mention that the differences 
observed in the position of the single particle levels in Fig. 
[3] has little impact on the quantum numbers of neighbor- 
ing odd-A nuclei as in the present mean field framework 
those quantum numbers have to be obtained after a self- 
consistent blocking mean field procedure and it is not 
enough to block the single particle orbitals of Fig. [3] as 
would be the case with a description based on a Nilsson 
diagram. Work to implement such blocking mechanism 
in the BCP case is under way and will be reported in the 
near future. 

In Fig. 5] we show the collective inertia B(Q 3 ) as- 
sociated with the octupole degree of freedom (see Eqs. 
(fTD) and (TTTjO and playing a central role in the collec- 
tive Schrodinger Hamiltonian of the previous section. As 
a consequence of the presence in its definition of a de- 
nominator with powers of the two quasiparticle energies 
+ E„, the collective inertia is roughly speaking in- 
versely proportional to the amount of pairing correla- 
tions (the pairing gap to be more quantitative) and di- 
rectly proportional to the effective mass of the interac- 
tion. The lower pairing correlations present in BCP1 
are not able to compensate for the higher effective mass 
and as a consequence the BCP1 inertias are higher than 
the Gogny D1S ones. Thus the energies obtained as a 
solution of the one dimensional collective Hamiltonian, 
which are roughly speaking proportional to the inverse 
of the square root of the collective mass (remember the 
standard harmonic oscillator formula relating the oscilla- 
tor's frequency oj with the spring constant and the mass 
to = y/k/m), are expected to reach lower values for BCP1 
than for D1S. It is also worth noticing that the peaks ob- 
served in the B(Qs) plots are related to regions of low 
pairing correlations as is easily deduced by comparing 
Fig. Hwith Fig. H 

In Fig [5] the zero point energy correction e(Qs) of Eq. 
|9] is given for the isotopes of radium considered and the 
BCP1 functional and the Gogny D1S force. The values 
°f ^{Qi) are correlated with the inverse of the collective 
inertia B(Qs), as can easily be noticed by comparing Figs 
@]aiid[SJ The range of variation is typically of around half 
an MeV in the interval of interest between Q3 — Ob 3 / 2 
and Q3 * 5b 3 / 2 and most of the nuclei considered al- 
though there are exceptions like the nucleus 216 Ra. The 
effect of the ZPE is to increase the depth of the octupole 
well for the lighter nuclei 216 ~ 220 Ra whereas it is the op- 
posite in all the heavier isotopes. The impact of this 
effect on the properties of the solutions of the CSE is not 
as pronounced as it could be thought because of the ef- 
fect of the collective masses (correlated to the behavior 
of the ZPE) that tends to cancel out the one of the ZPE. 
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Figure 4: The octupole collective inertia parameter B(Qs) en- 
tering the one dimensional collective Schrodinger Hamiltonian 
(see text for details) is shown as a function of the octupole mo- 
ment Q3 (in units of b 3 ^ 2 = 10 3 fm 3 ) for the radium isotopes 
considered. Full lines are for the calculations with the BCP1 
functional and dashed ones for the results obtained with the 
Gogny D1S force. 
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Figure 5: The octupole zero point energy correction e{Qz) of 
Eq. LS] is plotted as a function of the octupole moment Q3 (in 
units of b 3 ^ 2 = 10 3 fm 3 ) for the isotopes of radium considered. 
Full lines are for the calculations with the BCP1 functional 
and dashed ones for the results obtained with the Gogny D1S 
force. 
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Figure 6: The relevant quantities entering the Collective 
Schrodinger equation for the two cases considered, namely 
the calculation with the Gogny D1S force (right hand side 
panels) and the calculation with the BCP1 energy density 
functional (left hand side panels). In the upper panels, the 
octupole collective mass B{Qz) (already shown in Fig. Q is 
depicted as a function of the octupole moment Q3 (in units of 
^3/2 _ ^q3 £ m 3^ j n i ower p ane i S] the thick lines represent 
the HFB energy (dashed line) and the CSE potential energy 
(HFB energy minus zero point energy, full line). The thin 
lines correspond to the square of the collective wave functions 
corresponding to the lowest positive parity state \4> + (Q3) 2 
(full line) and lowest energy negative parity state \<j> - (Q3) | 2 
(dashed line). 



With the potential energy curve, the collective mass 
and the zero point energy correction, all the ingredients 
needed to solve the CSE are at our disposal. In Fig. [S] 
we have shown all those ingredients together in two plots 
corresponding to the results with the Gogny D1S force 
and BCP1 functional. In each of the plots we have de- 
picted in the lower panel the HFB energy curve (dashed 
line) as a function of the octupole moment Q3 and shifted 
to put the minimum at zero energy. The full curve closely 
following the dotted one is the potential energy entering 
the CSE that is obtained by subtracting the ZPE en- 
ergy correction to the HFB energy. As observed in the 
plot the collective potential energy is rather similar to 
the HFB energy. Also the HFB energies obtained with 
the D1S force and the BCP1 functional calculations are 
rather similar. In the same panel, the square of the col- 
lective amplitudes \<fi a (Q3)\ 2 for the lowest lying state of 
each parity are plotted. The negative parity amplitudes 
look rather similar in both calculations but this is not 
the case for the positive parity amplitude which is higher 
around Q 3 = for Gogny D1S than for BCPT The dif- 
ferent behavior of the positive parity amplitude is related 



to the different collective masses obtained in both calcu- 
lations and given in the upper panels. The Gogny D1S 
collective mass is much lower around Q 3 = than the one 
obtained with BCP1 and, as discussed in Ref. (Tol. Il7|. 
this enhances the collective amplitude around that value. 
As a consequence of the lower mass obtained with the 
D1S force the energy of the f "state computed after solv- 
ing the CSE is higher (200 keV) than the one obtained 
with the BCP1 functional (73 kcV). On the other hand, 
the effect on the B(E1, 1~ -> 0+) and B{E3,3~ ->■ 0+) 
transition probabilities is to yield smaller values for Df S 
than for BCP1 as the overlap between the positive and 
negative parity amplitudes is smaller in the later case. 
However, recalling the expression of the transition prob- 
abilities of Eqs. © and (J7J|, it is easy to realize the 
reduced impact of the region around Q3 = 0b 3 / 2 on the 
final quantities as each of the factors of the integrands, 
AKQs) and (Q3) prot (Q3), are zero for Q 3 = b 3 / 2 . 

The energies of the l - states and transition proba- 
bilities obtained by solving the CSE with the collective 
parameters deduced from the Gogny D1S and BCP1 cal- 
culations are depicted in Fig. [7] along with available ex- 
perimental values. For the energy of the 1~ states we 
observe that both the BCP1 functional and the Gogny 
Df S interaction reproduce quite nicely the experimental 
isotopic trend (see 0, IHJ and references therein) with a 
minimum around A=224. The very good reproduction 
of the experimental data in the calculation with Gogny 
D1S can be considered as accidental in the sense that the 
absolute values of the 1~ excitation energies crucially de- 
pend on the amount of pairing correlations (through the 
collective mass) which are not so well characterized at the 
mean field level. A more robust indicator of the quality 
of the results is the reproduction of the isotopic trend 
which is very good in both, D1S and BCP1, calculations. 
Concerning the B(E1, 1" — > + ) transition probabilities 
we observe a pronounced minimum around A — 224 in 
the two calculations that is a direct consequence of the 
behavior of the dipole moment as a function of the oc- 
tupole moment for different isotopes. This dip in the 
B(E1A~ — »■ + ) values is also observed experimentally 
(see jj, H3 and references therein) and is well reproduced 
by the Gogny Df S force and reasonably well by the BCP1 
functional. On the other hand, BCPf nicely reproduce 
the B(E1) of 226 Ra, whereas the Gogny Df S force result 
yields a too high value. Concerning the B(E3, 3~ — > + ) 
we observe a maximum around A = 224 which is corre- 
lated to the minimum in the energies of the 1 _ states. 
Both calculations re pro duce quite well the only experi- 
mental value known |33l . [34| . 

In order to get a more detailed understanding of the 
isotopic behavior of the B(E1 , f ~ — > + ), it is conve- 
nient to look at the behavior of the dipole moment Dq as 
a function of Q3 for the different isotopes considered. Ac- 
cording to Eq. [6l the value of the B(E1) transition prob- 
ability is proportional to the square of the average of the 
dipole moment over the whole (^interval and weighted 
with the product of the ground state positive parity col- 
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Figure 7: In the left hand side panel, the energies of the 1~ states, obtained by solving the ID collective Schrodinger Hamiltonian, 
are displayed as a function of the mass number for the radium isotopes considered. Theoretical results (lines, full for the 
results obtained with the BCP1 functional and dotted for the results obtained with the D1S force) are plotted along with the 
experimental data (boxes). In the middle panel the B(E1, 1~ — > + ) transition probabilities in W.u. are given as a function of 
the mass number for the radium isotopes considered. Finally, in the right hand side panel the B(E3, 3~ — > 0) in W.u. is given 
for the different radium isotopes considered. 
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Figure 8: Dipole moments Do (fin) are represented as a func- 
tion of the octupole moment Q3 (in units of b 3 ^ 2 = 10 3 fm 3 ) 
for the nine isotopes of Ra considered. Dashed lines cor- 
respond to the results obtained with the Gogny D1S force 
whereas the full line ones are for the results obtained with 
the BCP1 functional. 



therefore the dipole moment is always zero in that case. 
We observe that at the beginning of the isotopic chain the 
dipole moment increases monotonically as a function of 
Q^but its slope decreases with increasing neutron num- 
ber. For 224 Ra and also 226 Ra the slope is almost zero 
in the region from Q 3 — b 3 / 2 and up to Q3 = 5 b 3 / 2 , 
which is the region of interest where the collective wave 
function weight is different from zero. As a consequence, 
it is expected that the B(E1) has to reach a minimum for 
one of these isotopes. For 228 Ra and heavier isotopes the 
dipole moment in the region of interest decreases mono- 
tonically with a somehow constant slope what explains 
the increase in B(E1) in those isotopes as compared to 
224 Ra and 226 Ra as well as their almost constant value as 
a function of neutron number. The behavior of Dq with 
neutron number can be easily understood by looking at 
Fig. [3]whcre the SPE are plotted. There we observe how 
increasing the number of neutrons leads to the occupancy 
of more levels belonging to the high-j orbitals in/2 an d 
.715/2 ■ Due to the high total angular momentum value 
j of those orbitals, the spatial distribution of probabil- 
ity must have regions of large curvature that result in 
large values of (z) for those orbitals. Thus, increasing 
the number of neutrons increases the number of particles 
in those orbitals and the value of (z) neu t also increases, 
producing a decrease of Dq that is clearly seen in Fig. [5] 



lective wave function times the lowest negative parity 
one. The dipole moment entering Eq. [6] is represented 
as a function of Q3 in Fig. [8] for the Ra isotopes studied. 
Due to the good parity of the Q3 = b 3 / 2 solution the 
center of mass is located at the origin of coordinates and 



B. Neutron-rich barium isotopes 

Neutron rich barium isotopes (Z=56) with mass num- 
bers 142, 144 and 146 show several of the characteristics 
of octupole deformation in their ground states and yrast 
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Figure 9: The HFB mean field energies of barium isotopes 
are plotted as a function of the octupole moment Qg (in units 
of 6 3/2 = 10 3 fm 3 ) for the six even-even isotopes of barium 
(Z=56) considered. 



bands. Experiments [38] using the fragment yield of the 
252 Cf fission decay provided information on the yrast and 
negative parity rotational bands in these nuclei show- 
ing the typical alternating parity rotational band pat- 
tern representative of octupole deformed nuclei. For this 
reason we have performed calculations for the even-even 
barium isotopes with atomic numbers from A=140 up to 
A=150 to check the predictions of the BCP1 functional 
concerning octupolarity. Previous calculations with the 
Gogny D1S force in this region either at zero spin fl7l.l35l| 
or at high spins using the standard HFB cranking model 
[36| and even including parity projection (37j have been 
performed. In all the cases the agreement with experi- 
ment was satisfactory. 

In Fig. [9] we show for the Gogny D1S force (dotted line) 
and BCP1 functional (full line) the PEC corresponding 
to the six barium isotopes considered. Whereas in the 
Gogny D1S predictions it turns out that three isotopes 
have an octupole deformed minimum (namely, 144 Ba, 
146 Ba and 148 Ba) this is not the case for the BCP1 re- 
sults. However, in those nuclei the PEC calculated with 
the BCP1 functional are very flat around the Q3 = b 3 / 2 
minimum, which is a clear signature of a strong instabil- 
ity in the octupole degree of freedom. In addition, the 
depth of the octupole minima computed with the Gogny 
D1S force never exceed the 0.7 MeV found in the case 
of 144 Ba, which is a quite small height as compared to 
the typical energies of the vibrational octupole states. 
Therefore, the existence of the octupole minima can not 
be considered as conclusive. For the nuclei 142 Ba and 
150 Ba the results of both kinds of calculations show very 
flat curves around the Q3 = b 3 / 2 minimum indicating 
some degree of instability against the octupole degree of 



Figure 10: The particle-particle (pp) correlation energy de- 
fined as E pp — iTrAft; is plotted as a function of the octupole 
moment Q3 (in units of b 3 ^ 2 = 10 3 fm 3 ) for the barium iso- 
topes considered. Dotted lines correspond to the pp corre- 
lation energies of neutrons whereas full lines are meant to 
represent the same quantity for protons. Thick lines are used 
to depict the results of the calculations with the BCP1 func- 
tional whereas thin lines are for the results with the Gogny 
D1S force. 



freedom. Finally, the nucleus 1 °Ba is found to be rather 
stiff against octupole deformation in the two cases. 

In Fig [10] the particle-particle correlation energies are 
plotted as a function of the octupole moment. As in the 
case of the radium isotopes, we observe that the gen- 
eral trend of E pp both for protons and neutrons and the 
two interaction/functional considered is to decrease for 
increasing octupole moments in the relevant interval be- 
tween Q3 = b 3 / 2 and Q3 ~ 3.5b 3 / 2 ; a tendency that is 
also observed at higher values of the octupole moment in 
most of the barium isotopes studied. Based on the results 
of Fig. [Tolas well as the ones of Fig [2] for the radium iso- 
topes it is possible to say, in a very broad sense, that the 
onset of octupole deformation tends to quench pairing 
correlations and this quenching is bigger for larger values 
of the octupole moment. It is also noticed that in most 
of the cases the pairing correlation energies obtained in 
the BCP1 calculation are smaller than the ones obtained 
with the Gogny D1S force. 

In Fig. [TJJthe SPE spectrum of 144 Ba obtained with 
the two kinds of interactions/functionals is shown. The 
meaning of the different panels is the same as in Fig. [3] 
for 22 Ra. Here, the same comments made in the analysis 
and discussion of Fig. [3] for 224 Ra and regarding the Nils- 
son quantum number contents of the single particle wave 
functions are also in order. Coming back to the plot, 
we observe in the proton spectra the presence near the 
Fermi level of an occupied positive parity d 5 / 2 spherical 
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Figure 11: (Color online) Single particle energies (SPE) for the nucleus 144 Ba and computed with the BCP1 energy density 
functional (lower panels) and the Gogny D1S (upper panels) are plotted first as a function of the quadrupole moment Q2 (in b) 
up to the value corresponding to the self consistent minimum at Q3 = (Q2 = 3.37 b for the results obtained with the BCP1 
functional and 4.12 b for the ones with the Gogny D1S force) and from there one as a function of the octupole moment Qa(m 
b 3 / 2 ) (see the caption of Fig [3] for further details). 



orbital and a nearby empty negative parity hn/2 orbital. 
As discussed in the introduction, this is a characteristic 
property of the nuclear SPE for the octupole deformation 
to take place. In the SPE for neutrons we observe nega- 
tive parity orbitals coming from the fa / 2 sub-shell below 
the Fermi level and positive parity orbitals coming from 
the «i3/2 above, which is again a characteristic signature 
of octupole deformation. We also observe that the Fermi 
level both for protons and neutrons is located in a re- 
gion of low level density in the case of the Gogny force 
calculations, which is a required condition for octupole 
deformation to take place (Jahn- Teller effect). In the 
case of the calculations with BCP1 we observe that the 
Fermi level of neutrons is located in a region of moderate 
density of levels that could explain the lack of octupole 
deformed minima in the BCP1 calculations. 

As in the 224 Ra case, the position of the single particle 
levels cannot be used to assign the quantum numbers 



of neighboring odd-A nuclei in a fashion similar to the 
Nilsson diagram case. 

In Fig. [T^lthe energies of the 1~ states and the tran- 
sition probabilities obtained after solving the CSE are 
shown as a function of the mass number of the barium 
isotopes considered. In the left hand side panel the ex- 
citation energy of the 1~ state is plotted along with the 
available experimental data |38l - t4lT |. We observe how the 
isotopic trend is reasonably well reproduced although the 
absolute values of the energies are typically a factor of 
two larger than the experimental values. The BCP1 ex- 
citation energies are closer to the ones from Gogny D1S 
than in the radium isotopes calculations. As discussed 
previously, the values of the energies strongly depend 
upon the values of the collective mass in the vicinity 
of the minimum as well as on the height and width of 
the barrier separating the two minima with opposite oc- 
tupole deformation. Consequently, those excitation en- 
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Figure 12: Same as Fig. [7|but for the even-even barium isotopes with A in the range from 140 till 150. 



ergies depend on the interaction/functional used as well 
as the level of detail of the theoretical description (inclu- 
sion of pairing correlations, restoration of symmetries) 
and therefore the discrepancy with the experiment is not 
very relevant. Concerning the B(E1) values, we observe 
a dip in 146 Ba which is caused by the same effect as 
the dip in 224 Ra, namely the peculiar behavior of the 
dipole moment Dq{Q^) with the octupole moment and 
mass number. The BCP1 values for the B{E1) transition 
probabilities are systematically smaller than the Gogny 
D1S ones by almost one order of magnitude. As in the 
case of the excitation energies of the l - states, it has to 
be stressed that the isotopic trend is consistent in the 
two sets of calculations and both nicely reproduced the 
scarce experimental data (32l . [38| . Finally, in the right 
hand side panel of Fig. [12] the B(E3) transition proba- 
bilities are plotted. In this case and to our knowledge, no 
experimental data is available. The B(E3) values show 
a not very smooth behavior that, on the other hand, is 
inversely correlated with the excitation energies of the 1~ 
states. Therefore the bigger B(E3) values are obtained 
for the nuclei with the lower l - excitation energies. 



C. Fission valley properties of "Pu 

In this section the fission properties of 240 Pu, regard- 
ing the octupole contents of the mean field configura- 
tion in its way out to scission, are analyzed for the two 
BCP functionals [(| and the Gogny D1S force. This iso- 
tope has been chosen as a paradigmatic example of fis- 
sioning nucleus that has been thoroughly studied with 
the Gogny force [l4|. As it is customary, the theoreti- 
cal description of fission is based on the analysis of the 
potential energy curve (PEC) obtained by performing 
constrained mean field calculations with the quadrupolc 



moment Qi = (Q20) as constraining quantity. Details 
on the procedures involved can be found in the litera- 
ture Qj, EH, E3- The PECs obtained in this case for 
the nucleus 240 Pu are shown in the lower panel of Fig. 
1131 The range of quadrupole moments considered starts 
at sphericity (Q2 = 0) and goes up to values corre- 
sponding to very elongated configurations closed to fis- 
sion (Q2 = 150 b). At Q2 « 70 b the matter distribu- 
tion starts to resemble the form of two fission fragments 
connected by a neck as can be observed in the shapes 
given in the figure corresponding to the contour lines of 
the matter distribution at half density (see caption for 
further details). As the quadrupole moment increases, 
the distance between the fragments also increases and 
the width of the neck decreases. As a consequence, the 
energy roughly corresponds to the dominant Coulomb re- 
pulsion between the two incipient fragments. We observe 
how the results obtained with BCP and Gogny D1S are 
fairly similar with the position of the ground state min- 
imum and the fission isomer lying at roughly the same 
Q2 values. The potential energy curves obtained with 
the BCP1 and BCP2 functionals are very similar and 
are hardly distinguishable in Fig. Q2J The fission isomer 
obtained with Gogny D1S lies at an excitation energy al- 
most 2 MeV higher than the one obtained with the BCP 
functionals. It is also observed that the first and second 
fission barrier heights obtained with the Gogny D1S force 
are higher than the ones obtained from the BCP func- 
tionals. This is very likely a consequence of the surface 
coefficient in semi-infinite nuclear matter that is higher 
in Gogny D1S than in the BCP functionals (see the val- 
ues quoted in [Tl|). At this point and taking into account 
the differences in the PECs it can be concluded that the 
predictions for the spontaneous fission half lives obtained 
with the BCP functionals and Gogny D1S force are going 
to be higher for Gogny D1S than for BCP. However, a 
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Figure 13: (Color online) In the lower panel the axially sym- 
metric fission path (mean field energy versus axial quadrupole 
moment in barns) for the nucleus 240 Pu and the three inter- 
actions/functionals considered (BCP1 full line, BCP2 dashed 
line, Gogny D1S dotted line). The shapes depicted corre- 
spond to the half-density contour line of the actual matter 
density distribution and are plotted such that their symme- 
try axis is located at the corresponding Q20 value. In the 
upper panel the octupole and hexadecapole moments of the 
corresponding mean field states. Apparently, there are only 
two curves but this is because the results for the three func- 
tionals/interactions lie one on top of each other. 

definitive answer to this question can not be given un- 
til the effect of triaxiality has been incorporated into the 
calculations because it is well known that triaxiality can 
have a strong impact in the first barrier height. Also, it 
has to be kept in mind that the collective mass along the 
Q2 collective degree of freedom, and entering the WKB 
formula used to estimate fission half lives, can be sub- 
stantially different when computed with the BCP func- 
tionals or the Gogny D1S force. Therefore, the detailed 
discussion of the fission half lives obtained with the BCP 
functionals is deferred to a more detailed study of fission 
properties obtained with this class of functionals. In this 
paper, devoted to octupole deformation, it is enough to 



confirm that the shapes of the nucleus in its way down to 
fission are essentially the same irrespective of the func- 
tional/interaction used as can be seen in the upper panel 
of Fig. [T31 In this plot, the octupole and hexadecapole 
moments are depicted as a function of Q2 and, as can be 
observed, the curves for different interactions/ functionals 
are indistinguishable of each other. This fact implies that 
the mass distribution close to scission (as depicted in Fig. 
[TBI through the half density contours) is the same irre- 
spective of the interaction/functional used in the calcula- 
tion and therefore the predictions of the fission fragment 
mass distributions obtained at the mean field level with 
Gogny D1S force and BCP functionals should coincide. 



IV. CONCLUSIONS 

We have explored the octupole degree of freedom in 
two sets of isotopes with the newly postulated BCP func- 
tionals. The agreement found with both experiment and 
the benchmark results obtained in the same framework 
with the Gogny D1S interaction, gives us confidence on 
the good properties of the BCP functionals concerning 
odd parity multipole moments. In addition, the mat- 
ter distribution of the fissioning nucleus 240 Pu, which 
strongly depend upon the response of the system to oc- 
tupole perturbations, is found to be essentially the same 
in the three calculations performed implying thereby that 
the BCP functionals and the Gogny D1S force are equally 
well suited in that respect. Taking into account the mi- 
croscopic origin of the BCP functionals it is comforting 
and encouraging to observe its good performance in prop- 
erties like octupolarity that belong to the realm of finite 
nuclei. 
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